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We review the present status of the Baikal Neutrino Project and present results on upward going atmospheric 
neutrinos, results of a search for high energy extraterrestrial neutrinos as well as preliminary results of searching 
for acoustic signals from EAS in water. We describe the moderate upgrade of NT-200 planned for the next years 
and discuss a possible detector on the Gigaton scale. 


1. DETECTOR AND SITE 

The Baikal Neutrino Telescope is operated in 
Lake Baikal, Siberia, at a depth of 1.1 km. The 
present stage of the telescope, NT-200 [I|, was 
put into operation at April 6th, 1998 and consists 
of 192 optical modules (OMs). An umbrella-like 
frame carries 8 strings, each with 24 pairwise ar¬ 
ranged OMs. Three underwater electrical cables 
and one optical cable connect the detector with 
the shore station. 

The OMs are grouped in pairs along the strings. 
They contain 37-crn diameter QUASAR - photo 
multipliers (PMs) which have been developed spe¬ 
cially for our project [||. The two PMs of a pair 
are switched in coincidence in order to suppress 
background from bioluminescence and PM noise. 


A pair defines a channel. 

A muon-trigger is formed by the requirement of 
> N hits (with hit referring to a channel) within 
500 ns. N is typically set to 3 or 4. For such 
events, amplitude and time of all fired chan¬ 
nels are digitized and sent to shore. A separate 
monopole trigger system searches for clusters of 
sequential hits in individual channels which are 
characteristic for the passage of slowly moving, 
bright objects like GUT monopoles. 

Lake Baikal deep water is characterized by an 
absorption length of L a bs (480 nm)=20 =24 m, a 
scattering length of L s =30 =70 m and a strongly 
anisotropic scattering function f(0) with a mean 
cosine of the scattering angle cos(0) = 0.85 = 0.9. 
Fig.l shows the cascade detection volume and 
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Figure 1. Detection volume (left) and detec¬ 
tion area (right) of a single BAIKAL OM for neu¬ 
trino induced high energy cascades and high en¬ 
ergy muons, respectively. 


the muon detection area of a single BAIKAL 
OM. Here, we define detection area and detec¬ 
tion volume by the condition that the mean num¬ 
ber of photoelectrons has to be >1. In contrast 
to underground detectors, open configurations in 
highly transparent media like water or ice allow 
to observe a huge volume beyond their geometri¬ 
cal boundaries. The detection volume of an OM 
rises from 1-10 5 m 3 for 1 TeV to 7-10 7 m 3 for 1 
EeV cascade energy. 

Here we present results on upward going at¬ 
mospheric neutrinos, a search for high energy ex¬ 
traterrestrial neutrinos, and preliminary results 
of a search for acoustic signals from EAS in water. 
Also an upgrade of NT-200 towards the Gigaton 
scale detector is dicussed. 

2. A SEARCH FOR UPWARD GOING 
ATMOSPHERIC NEUTRINOS 

The signature of neutrino induced events is a 
muon crossing the detector from below. The re¬ 
construction algorithm is based on the assump¬ 
tion that the light radiated by the muons is emit¬ 
ted under the Cherenkov angle with respect to 


the muon path. We don’t take into account light 
scattering because the characteristic distances for 
muons induced by atmospheric neutrinos which 
trigger the detector do not exceed 1-^2 scattering 
lengths of light in Baikal water (see above). 

The algorithm uses a single muon model to re¬ 
construct events. We apply procedures to reject 
hits, which are very likely due to dark current 
or water luminosity as well as hits which are due 
to showers and have large time delays with re¬ 
spect to expected hit times from the single muon 
Cherenkov light. 

Determination of the muon trajectory is based 
on the minimization of a y 2 function with respect 
to measured and calculated times of hit channels. 
As a result of the \ 2 minimization we obtain the 
track parameters (0, (f> and spatial coordinates). 
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Figure 2. Experimental angular distribution of 
events reconstructed and accepted as upward go¬ 
ing muons (solid crosses). The MC expectation 
for atmospheric neutrinos is shown for the case 
of no oscillation (solid line) and for oscillations 
with maximal mixing and Am 2 = 3- 10 -3 eV 2 
(dashed line). (MC distribution do not include 
the background due to atmospheric muons.) 

Since the algorithm yields a large fraction of 
events which are reconstructed as upward going 
muons, wrongly reconstructed tracks have been 
removed by a series of quality cuts j|||]. 

The efficiency of reconstruction and “fake” 
event rejection and the correctness of the MC 
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background estimation have been tested with a 
sample of 1.1-10 8 MC-generated events from at¬ 
mospheric muons (that is twice the whole ex¬ 
perimental sample, see below) and with MC- 
generated events due to atmospheric neutrinos. 
Muons passage through rock and water has been 
generated by transport code MUM ||. 

For the MC background sample we found 3 
surviving events reconstructed as upward going 
muons and passing the quality cuts. The anal¬ 
ysis of these events shows that they are com¬ 
plex and contain one or two muons with show¬ 
ers (shower energy > 100 GeV). These events are 
reconstructed close to horizon (cos(0) > —0.3). 

Data taken with NT-200 between April 1998 
and February 1999 cover 234 days life time. From 
the total of 1.67-10 8 triggers (IV^ >4) recorded 
in 234 days, 5.3-10 7 events with Aqq^ >6 at > 3 
strings have been selected for this analysis. These 
events have been reduced to 34 upward tracks. 
The MC estimations of the number of expected 
upward going events due to neutrino induced 
muons give 34 and 29.4 events in case of absence 
of neutrino oscillations and including oscillations 
with maximal mixing angle and Am 2 =3-10~ 3 
eV , respectively. Angular distributions of both 
experimental and MC samples are presented in 

Fig# 

Fig. H shows the celestial distribution of selected 
upward going muons with NT-200 in galactic co¬ 
ordinates. 
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Figure 3. The celestial distribution of the exper¬ 
imental neutrino sample in galactic coordinates. 
The thick solid line indicates the horizon in equa¬ 
torial coordinates. 


3. A SEARCH FOR EXTRATERRES¬ 
TRIAL HIGH ENERGY NEUTRINOS 

The used search strategy for high energy neu¬ 
trinos relies on the detection of the Cherenkov 
light emitted by the electro-magnetic and (or) 
hadronic particle cascades and high energy muons 
produced at the neutrino interaction vertex in a 
large volume around the neutrino telescope HQ. 
A cut is applied which accepts only time patterns 
corresponding to upward traveling light signals. 

For this analysis we used 1.56T0 5 events with 

A hit >10. 



Figure 4. Distribution of hit channel multiplicity; 
dots - experiment, hatched boxes - expectation 
from brems and hadronic showers produced by 
atmospheric muons. 


Fig.4 shows the A^n distribution for ex¬ 
periment (dots) as well as the one expected 
for the background from brems- and hadronic 
high energy showers produced by atmospheric 
muons (boxes). The experimental distribution is 
consistent with the background expectation for 
ATjjif >18. The discrepancy at lower A^p. val¬ 
ues is caused by the contribution of atmospheric 
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muons close to horizon as well as by low energy 
showers from e^e - pair production. These ef¬ 
fects have not been included in the present MC 
sample, but have been taken into account in sepa¬ 
rate simulations which are consistent with exper¬ 
imental data. No statistically significant excess 
over background expectation from atmospheric 
muon induced showers has been observed. Since 
no events with iV^jf >45 are found in our data, 
we can derive upper limits on the flux of high en¬ 
ergy neutrinos which would produce events with 
N hk >50. 

The detection volume V e jf for neutrino pro¬ 
duced events with TVjjU >50 which fulfill all trig¬ 
ger conditions was calculated as a function of neu¬ 
trino energy and zenith angle 9. V e ff rises from 
2-10 5 m 3 for 10 TeV up to 6T0 6 m 3 for 10 4 TeV 
and significantly exceeds the geometrical volume 
V g « 10 5 m 3 of NT-200. 

Given an E~ 2 behaviour of the neutrino spec¬ 
trum and a flavor ratio (v e -\-v e ) '■ (pu+^v) = 1 : 2, 
the combined 90% C.L. upper limit obtained with 
the Baikal neutrino telescopes NT-200 (234 days) 
and NT-96 g (70 days) is: 

$ (l/e+I ; e) E 2 < (1.3>1.9)T0 _6 cm _2 s _1 sr _1 GeV(l) 

where the upper value allows for the highest light 
scattering observed over many seasons. 

Fig.5 shows the upper limits on the isotropic 
diffuse neutrino flux obtained by BAIKAL (this 
work), AMANDA [g, MACRO [g and FRE- 
JUS (lO) (triangle) as well as the atmospheric 
conventional neutrino fluxes ]Il| from horizontal 
and vertical directions (upper and lower curves, 
respectively) and the atmospheric prompt neu¬ 
trino flux |l2j (curve labeled v pr ). Also shown is 
the model-independent upper limit on the diffuse 
high energy neutrino flux obtained by Berezin¬ 
sky jU| (curve labeled ’B’), and predictions for 
diffuse neutrino fluxes from Stecker and Salamon 
fj (’SS’) and Protheroe |li| (’P’). Curves la¬ 
beled ’MPR’ and ’WB’ show the upper bounds 
obtained by Mannheim et al. [|T(| as well as the 
upper bound obtained by Waxman and Bahcall 
fl7|| , respectively. Curves labeled ’M(GRB)’ and 
’WB(GRB)’ present the upper bounds for diffuse 
neutrino flux from GRBs derived by Mannheim 
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Figure 5. Experimental upper limits on the neu¬ 
trino fluxes as well as flux predictions in different 
models of neutrino sources (see text). 


fni and Waxman and Bahcall jP)]]. The curve la¬ 
beled ’TD’ shows the prediction for neutrino flux 
from topological defects due to specific top-down 
scenario BHS1 [pof . 

Our combined 90% C.L. limit at the W - reso¬ 
nance energy is: 

< (1.4 >1.9) x 10 _19 cm _2 s _1 sr _1 GeV _1 (2) 

dE p 

and is given by the rectangle in Fig.5. 

4. SEARCH FOR ACOUSTIC SIGNAL 
FROM EAS CORE IN WATER 

In 2002 the Baikal collaboration has reported 
the observation of characteristic bipolar acoustic 
signal in time coincidence with an EAS detected 
by a scintillator array g. The scintillator array 
was placed on ice cover and the hydrophone in 
water, at the depth of 5 m and at a distance of 90 
m from the array center. In the two years since 
then experiments with acoustic hydrophones trig¬ 
gered by the EAS array have been continued. 
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In both experiments the EAS array consisted of 
seven scintillator counters. Six counters formed a 
circle with radius 80 m in 2001 (50 m in 2002). 
The seventh counter was placed in the center of 
the circle. Measurements of the number of rel¬ 
ativistic particles in each counter and the rela¬ 
tive time delays allow reconstructing the shower 
parameters: arrival direction, core position and 
total number of electrons. 

In the experiment of March/April 2001, four 
hydrophones were located at a distance of 34 m 
from the center of the EAS array and formed a 
square. Ten runs of measurements have been car¬ 
ried out. Hydrophones were placed in water at a 
depth of 4 m in seven runs, while they were lo¬ 
cated in air, on the ice cover, in three runs. 

In the second experiment (March/April 2002), 
a more compact acoustic antenna was used. 
Three hydrophones were placed at the corners of 
an equilateral triangle at a distance of 1 m each 
from other. All three hydrophones were in water 
at a depth of 4 m. The fourth hydrophone was 
located below the triangle center at a depth of 8 
m. 

On top of our own anntenas, hydrophones from 
an ITEP/Moscow group have been operated in 
both experiments. ITEP’s antennas had four hy¬ 
drophones of different type, placed on a vertical 
string and separated from each other by 5 m [plf . 

In order to record an acoustic signal, a multi 
channel digitizer was used. The EAS array trig¬ 
gered the data acquisition system and started the 
digitizers. The digitizers recorded the signal from 
each hydrophone during 0.1 sec with a step width 
of 2 /isec. So, for each shower there are four time 
series consisting of 50 000 measurements of the 
amplitude of the acoustic signal. 

The analysis of the time series shows that apart 
from Gaussian noise, which can be parameterized 
by average amplitude and standard deviation a, 
there are large fluctuations lasting from tens up to 
hundreds microseconds. These ’’perturbations” 
could be due to real acoustic signals produced by 
some source. Selecting ’’perturbations” with a 
duration between 20 and 100 /isec and an aver¬ 
age amplitude more than 2<j, we have found in 
average 50 ’’perturbations” in each channel per 
one event. Probably, they are mainly due to ice 


cracking. 

Combining ’’perturbations” recorded by the 
four hydrophones and using the known sound ve¬ 
locity in Baikal water, 1405 m/sec, one can re¬ 
construct the position of the source in space and 
time. We have found 46 events within a time 
window +/-0.5 msec relative to the EAS trigger, 
when the hydrophones were placed in water, and 
13 events in the case of hydrophone location in 
air. 

The expected background was estimated as¬ 
suming that the sound sources are independent 
and uniformly distributed in ice. Simple Monte 
Carlo calculations predict 26 and 12 events re¬ 
spectively. One can say there is excess of events 
within time window +/- 0.5 msec at the level 
of 2cr in the case of hydrophones located in wa¬ 
ter. The analysis of the 2002 experiment is in 
progress. 

5. NT-200+ AND BEYOND 

Recently derived upper limits on 
(AMANDA) and v e (BAIKAL) fluxes are about 
E 2 F(y) k, 10 -6 cm -2 s -1 sr -1 GeV and be¬ 
long to the region of most optimistic theoretical 
predictions. A flux sensitivity at the level of 
E 2 F{y ) < 10~ 7 cm -2 s -1 sr -1 GeV which would 
test a variety of other models, requires detection 
volumes of a few dozen Mtons. 

We envisage an upgrade of NT-200 to the scale 
of a few dozen Mton by three sparsely instru¬ 
mented distant outer strings. The basic princi¬ 
ple will be the search for cascades produced in a 
large volume below NT-200. This configuration, 
christened NT-200+, will not only result in an 
increased detection volume for cascades, but also 
allow for a precise reconstruction of cascade ver¬ 
tex and energy within the volume spanned by the 
outer strings. 

A schematic view of NT-200+ is shown in Fig.6. 
It will comprise the neutrino telescope NT-200 
itself as well as three 140 nr long outer strings with 
3 pairs of OMs spaced vertically by 70 m, and two 
pairs of OMs spaced by 6 m and arranged at one 
of the eight strings of NT-200, 140 m below it. 
The outer strings are arranged at a distance of 100 
nr around NT-200 at the edges of an equilateral 
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Figure 6 . Sketch of NT-200+. 


triangle. Their top OMs are located at the level 
of the bottom OMs of NT-200. A water volume 
of 4.4 • 10 6 m 3 is surrounded by the outer strings 
and NT-200. 

Moderate event selection requirements allow 
to achieve a large detection volume for neutrino 
events and to suppress background effectively. 
The detection volumes for isotropic v e and 
fluxes are shown in Fig.7. The value of 14// for v e 
induced events rises from 7-10 5 m 3 for 100 TeV up 
to 4-10 7 m 3 for 10 6 TeV. Normalized energy dis¬ 
tributions of expected event rates from v e and v )J: 
fluxes following an inverse power law with spec¬ 
tral index 7 = 2 are presented in Fig. 8 . Most of 
the expected events would be produced by neu¬ 
trinos from the energy range 10 2 4-10 5 TeV, with 
a mean energy around 1 PeV. 

Assuming 7 = 2 and flavor ratio (v e + u e ) : 



lg(E,/TeV) 


Figure 7. Detection volume of NT-200+ for v e 
and v^ events which survive all cuts. 


(r/j + 14) = 1:2 and a v e flux of 

<f>(j, e+! 7 e ).E 2 = 3.5 • 10 _ 7 cm _ 2 s _ 1 sr _1 GeV (3) 

a one year expectation would yield 2.4 events. 

MC simulations show that the detection vol¬ 
ume of NT-200+ for PeV cascades varies only 
moderately, if NT-200 as the central part of NT- 
200-1- is replaced by a single string of OMs. Fig.9 
gives the detection volume for different config¬ 
urations as a function of cascade energy. The 
standard configuration of NT-200+ is marked by 
empty rectangles. The other configurations com¬ 
prises a single string instead of NT-200: a stan¬ 
dard string of 70 m length and 24 OMs (filled rect¬ 
angles), a half string with 12 OMs covering 35 m 
(dots), and a 70 m long string sparsely equipped 
with 12 OMs (triangles). The configuration with 
the long 12-OM string shows an energy behaviour 
very close to the one of NT-200+. For neutrino 
energies higher than 100 TeV such a configura¬ 
tion could be used as a basic subarray of a Giga- 

















7 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

lg(E„/TeV) 


<t> 


10 ' 


o 

> 

c 

_o 

o 

a) 

-t-- 

<t> 

~o 


10 ' 


10 - 


10 


- 

V, 




** 


: 


r 

..g 

it 

i 






rf 

E 

i=E 

cflS 1 

* 


- NT2 

, - long 

00 

strin 

T9 

.. ^ . 

OM 

) 

..... # 

1 “ 

• 

....A... 

■ 

▲ 

* 

▲ 



— long 

, _ half ! 

strin 

trin 

g(i: 

? (12 

OM 

OM 

) 

) 

:▲ 

mi 

mi 




mi 




1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

lg(E„/TeV) 


Figure 8. Energy distribution of expected events 
induced by diffuse v e and fluxes. 


ton Volume Detector (GVD). Rough estimations 
show that 0.7 4- 0.9 Gton detection volume for 
neutrino induced high energy cascades may be 
achieved with about 1300 OMs arranged at 91 
strings. A top view of GVD as well as sketch 
of one basic subarray are shown in Fig.10. The 
physical capabilities of GVD at very high energies 
cover the typical spectrum of cubic kilometer ar¬ 
rays. We are presently working on simulations to 
optimize the response for TeV muons, maintain¬ 
ing at the same time the cubic kilometer scale for 
cascades with energy above 100 TeV. 

6. CONCLUSIONS AND OUTLOOK 

The deep underwater neutrino telescope NT- 
200 in Lake Baikal is taking data since April 1998. 
Using the first 234 live days, 34 neutrino induced 
upward going muons have been selected. The lim¬ 
its on the diffuse high energy (i> e + D e ) flux as well 
as v e flux at the W-resonance energy have been 


Figure 9. Detection volumes of different configu¬ 
rations. 


derived. Also limits on an excess of the muon flux 
due to WIMP annihilation in the center of the 
Earth and on the flux of fast magnetic monopoles 
have been obtained. The additional analyzes of 
the experimental data accumulated during 1999- 
2002 will allow to decrease these limits by a factor 
3-4. 

Within the next few years we plan to upgrade 
the Baikal neutrino telescope to the 10 Mton 
detector NT-200+ with a sensitivity of approx¬ 
imately 3.5 • 10~ 7 cm _2 s _1 sr _1 GeV for a diffuse 
neutrino flux within the energy range 10 2 TeV 4- 
10 5 TeV. NT-200+ will search for neutrinos from 
AGNs, GRBs and other extraterrestrial sources, 
neutrinos from cosmic ray interactions in the 
Galaxy as well as high energy atmospheric muons 
with E^ > 10 TeV. 

In parallel to this short term goal, we started 
research & development activities towards a Gi- 
gaton Volume Detector in Lake Baikal. 

This work was supported by the Russian Min- 




























































Figure 10. Top view of GVD as well as sketch 
of one of its subarrays. 
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